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FLUTTER DESIGN OF STIFFENED-SKIN PANELS 
FOR HYPERSONIC AIRCRAFT 

By Herman L. Bohon, Melvin S. Anderson, 
and Walter L. Heard, Jr. 

Langley Research Center 

SUMMARY 

Factors which control the flutter design for exterior skin panels of hypersonic vehi- 
cles are reviewed. From a typical trajectory the most critical flight regime is identified 
in the supersonic range. Because past studies have shown that support flexibility can 
seriously reduce the flutter margin of panels such as those proposed for hypersonic vehi- 
cles, conservative design curves for such panels with flexible supports are presented. 
Means of greatly improving flutter margins by the addition of torsional restraint at the 
panel edges are also indicated. 


INTRODUCTION 

External skin panels of hypersonic aircraft will consist of various stiffened panels 
that may be subject to flutter in the supersonic and hypersonic operating range. Studies 
indicate that these panels will be of more complicated construction than current practice 
in order to withstand the hostile environment with minimum mass. Although a wealth of 
experimental data on the flutter of flat panels has been accumulated over the past few 
years (ref. 1), only a small portion of the data was obtained from stiff ened- skin panels. 
Further, only within the past 5 years have flutter data been generated at Mach numbers 
higher than Mach 5 (ref. 2), and nearly all these data have been obtained from isotropic 
panels. Nevertheless, data on isotropic panels suggest that panel flutter boundaries 
based on supersonic aerodynamic theory can be extended into the hypersonic range. Fur- 
ther, consideration of the flight envelope for proposed hypersonic cruise vehicles indicates 
the most critical portion of the flight from a flutter standpoint is in the supersonic- speed 
range. It is in this speed range that most flutter data on orthotropic panels have been 
obtained. Thus, the use of these data and recent theory derived for the supersonic- speed 
range permit the assessment of panel flutter design for hypersonic vehicles. In the pres- 
ent paper the available experimental and theoretical information is used to identify the 
critical flutter problems and a simplified approach for hypersonic vehicle panel design is 

given. 



SYMBOLS 


The units used for the physical quantities in this paper are given both in the U.S. 
Customary Units and in the International System of Units (SI). Factors relating the two 
systems are given in the appendix for the units used in the present investigation. 
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FLUTTER DESIGN CONSIDERATIONS FOR PANELS 
OF HYPERSONIC VEHICLES 


External skin panels of hypersonic vehicles will be exposed to a wide range of 
dynamic pressure, Mach number, and temperature. The impact of some of these param- 
eters on panel flutter design is discussed in the following sections. 

Effects of Mach Number 

Experimental flutter data generated from flat isotropic panels at high Mach num- 
bers (see ref. 2) have shown reasonable agreement with trends predicted by using small- 
deflection plate theory with quasi-static aerodynamics, as illustrated in figure 1, where 

the parameter q/D^M^ 2 - 1 is plotted as a function of local Mach number. The solid 
line is the flutter boundary, which separates the "flutter" region (above the line) from 
the "no flutter" region (below the line). This boundary, based on quasi-static aerodynam- 
ics, is constant with respect to Mach number and extends into the hypersonic range. The 
symbols represent experimental flutter data from reference 2, obtained from isotropic 
panels over the range of Mach numbers from 2 to 6. The agreement (considered good for 
panel flutter) of the experimental data with theory over the entire Mach number range in 
figure 1 suggests that the flutter boundary may be applied further in the hypersonic range 
than shown in figure 1 with reasonable success. The fact that the flutter boundary is not 
dependent on Mach number also suggests that panels on a hypersonic vehicle will be most 

susceptible to flutter wherever the flight path produces the maximum value of q/\/^ 2 - 1 
(if the panel bending stiffness D is assumed to be constant). Thus, trajectories for the 
hypersonic vehicle need to be examined to determine where the maximum value occurs. 
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Figure 2 shows the variation of the parameter - 1 with free-stream Mach num- 

ber for the wing lower surface of a hypersonic cruise vehicle in a typical ascent trajec- 
tory. In the parameter, the dynamic pressure and the Mach number are local conditions. 
The lower curve represents the local flow conditions over a flat surface at zero angle of 
attack, whereas the upper curve corresponds to flow conditions on the wing lower surface 
at angle of attack of 3° during the early part of the flight and to an angle of attack of 5° 
during the final climb at constant dynamic pressure. The point where a panel is most 

susceptible to flutter is at the maximum value of q/^M^ 2 - 1, which occurs on this flight 
path in the vicinity of Mach 3. Several other advanced vehicle trajectories, examined in 
reference 3, have led to the same conclusion — that the supersonic range is most critical 
from a flutter standpoint. 

The results of figure 2 would be modified slightly if the effect of elevated tempera- 
tures on the panel bending stiffness D is considered. However, since the design condi- 
tion for flutter is in the supersonic range where temperatures are generally low relative 
to material capability, the effect of variable D would be small. 

No mention has been made of the transonic range where, as is well known from 
experimental data, reductions in the critical value of f" occur. (See refs 4 

5, and 6.) Note, however, from figure 2 that the value of q/^M^ 2 - 1 encountered on 
the flight path also decreases in the transonic range; this reduction is generally greater 
than that obtained experimentally for flutter. Consequently, flutter in the transonic range 
is not pursued further in this paper. 

Effects of Panel Construction and Boundary Conditions 

Most experimental data available on stiffened panels suitable for a hypersonic vehi- 
cle have been obtained in the supersonic range where problem areas having special signif- 
icance with regard to panel flutter design have been well delineated. Since supersonic 
flutter analyses and experiments are pertinent to the hypersonic range, the effects of 
these problems on hypersonic panel design can be determined. 

Panel configurations appropriate for use on the hypersonic vehicle include waffle- 
grid, honeycomb core sandwich, or corrugation-stiffened panels with either chordwise or 
spanwise orientation. One of the primary design problems associated with such panels is 
the method of attachment to the supports. The sensitivity of panel flutter to support con- 
ditions is illustrated in figure 3, which shows two typical methods of attaching corrugation- 
stiffened panels at the ends of the corrugations. The panels which were identical except 
for support conditions were tested at Mach 3 with the corrugations alined normal to the 
airflow; thus, the weaker of the panel bending stiffnesses Dj is in the airflow direction. 
(Test details are reported in ref. 7.) In one test the corrugations at the edges were 
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attached to an angle clip; the panel fluttered at a dynamic pressure q of 3400 psf 
(163 kN/m2). The same panel, when tested with the clip removed (that is, with the cor- 
rugations unsupported), fluttered at a dynamic pressure q of 540 psf (25.8 kN/m2). 
Thus, the flutter margin was reduced by more than a factor of 6 simply by a relaxation 
of the stiffness of the supports at the ends of the corrugations. In order to determine 
analytically the flutter mechanisms involved in these and other similar tests, the problem 
has been formulated, in the next section, to account for arbitrary deflectional stiffness at 
the supports. 


ANALYTICAL RESULTS FOR ORTHOTROPIC PANELS 

An orthotropic panel and the coordinate system are shown in figure 4. The edges 
at x = 0 and x = a are simply supported. The edges at y = ±b/2 are supported by 
deflectional springs with a spring constant k D per unit length. These edge conditions 
may account for flexibility of the supports as well as for local deformations of the cross 
section at the support. The panel has been analyzed for airflow (at Mach number Mj) 
parallel either to the x-direction or the y-direction; the lateral loading due to the air 
forces is given by the two-dimensional static approximation. For details of the analyti- 
cal development, see reference 8. Numerical results obtained are applicable to any com- 
bination of bending and twisting stiffnesses and panel geometry, but the main part of the 
discussion herein applies to orthotropic panels for which Dj is much less than Di2> 
which is, in turn, less than D 2 * 

Spring Supports on Leading and Trailing Edges 

When the airflow is parallel to the y-direction (direction of strong flexural stiff- 
ness), the leading and trailing edges are supported by deflectional springs. Results, 
obtained from a closed-form solution for zero midplane stress, are shown in figure 5 on 
a plot of X cr , the dynamic pressure parameter at flutter, as a function of stiffness - 
geometry parameter -Ay for various values of nondimensional spring constant Kq. 

It should be noted that the panel dimension b is in the airflow direction. Flutter bound- 
aries are shown for finite values of spring constant Kq of 1, 5, and 10. The upper 
curve for Kq = °° is the flutter boundary for the panel with all edges simply supported. 
Values of Kq as low as 1 are not unreasonable, and the figure shows that support flexi- 
bility can reduce the flutter dynamic pressure q by as much as a factor of 10 from the 
simple- support value. At low values of -A y (panels with corrugations alined with the 
airstream), the curves become horizontal; hence, the flutter boundaries become indepen- 
dent of the panel width a and of the support conditions along the streamwise edges. In 
fact, the panel behaves as if it were infinitely long in the weak bending direction (that is, 
the x-direction). 
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Spring Supports on Streamwise Edges 

When the airflow is parallel to the x-direction (that is, normal to the direction of 
corrugations), the spring supports are at the streamwise edges. A closed-form solution 
could not be found for this support condition; instead, a Galerkin solution which used up 
to 50 natural panel vibration modes to assure convergence was obtained. The analysis is 
presented in reference 8, and some numerical flutter results are shown in figure 6. 

The parameters in figure 6 are slightly modified compared with those of figure 5 
as a result of the different direction of airflow. The cube root of \* T is plotted as a 
function of -A x for various values of K D . However, for this analysis the results 
are also a function of the parameter C, which is a ratio of panel stiffnesses shown in 
the figure. A representative value of C = 7 was used for these calculations. The 
curves for finite values of Kp are bounded above by the limiting curve for Kj) = «■ 

(that is, for simply supported edges), which is obtained from an exact solution, and below 
by the limiting curve for K D = 0 (that is, for free edges), which is also obtained from 
an exact solution. The cusps in the curves are caused by changes in critical modes which 
coalesce for flutter. The numerals in the regions separated by the dashed curves indicate 
the modes that coalesce to give the critical flutter boundary. 

For corrugation- stiffened panels, the ratio is very large; therefore, the 

region of primary interest is the far right-hand section of figure 6. The results show 
that flexible supports may cause large reductions in A* r . In this region, the flutter 
boundaries become horizontal and the flutter dynamic pressure q becomes independent 
of the panel length a; thus, flutter results are independent of details of the boundary con- 
ditions at the leading and trailing edges. 

Effect of Orientation of Corrugations 

The effect of orientation of the corrugations with respect to the airstream is shown 
in figure 7 where a direct comparison is made between the results for flow parallel to 
and the results for flow perpendicular to the orientation of the corrugations. The com- 
parison is made for a value of the ratio a/b equal 1. The dashed curves are from the 
exact analysis made with the corrugations oriented parallel to the airstream; the solid 
curves are from the modal analysis made with the corrugations oriented normal to the 
airstream. As noted previously, the modal results are functions of the stiffness param- 
eter C; curves are shown only for C = 7. For orthotropic panels £>12/^2 small, 
and the results show a pronounced adverse effect of weak supports on the streamwise 
edges. Additional effects of panel orientation can be obtained from reference 9 where 
flutter boundaries for Kq = « are presented for arbitrary angles of orientation and 
for several combinations of panel stiffnesses. 
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DESIGN BOUNDARIES 


As indicated in the discussion of figures 5 and 6, panels of interest for hypersonic 
vehicles generally fall in the regions where the flutter boundaries become horizontal. 

These regions of minimum values may provide good, and often conservative, estimates 
of the flutter parameter for orthotropic panels and should be very useful for design pur- 
poses. The plot of such minimum values of critical flutter parameter is shown in fig- 
ures 8 and 9 where A cr and A* r , respectively, are plotted as a function of K D over 
the range from K D = 0 (free edges) to Kq = «° (simply supported edges). A family of 
curves is shown in figure 9 for several values of the stiffness parameter C. The appli- 
cation of these curves by the design engineer requires a knowledge of the panel stiffnesses 
and geometry and a reliable estimate of the deflectional stiffness at the supports. Deflec- 
tional stiffness may be obtained experimentally from a vibration survey of the panel. 

To date, no experimental flutter data have been obtained on orthotropic panels ori- 
ented with the maximum flexural stiffness in the airstream direction to verify or dispute 
the trends shown in figure 8. On the other hand, all experimental flutter data obtained to 
date on orthotropic panels have been with the maximum flexural stiffness normal to the 
airstream, as represented by the theoretical model for the curves in figure 9. A corre- 
lation of data from several experimental investigations with the present theory has been 
made and the results are shown in figure 10. The solid curves are repeated from figure 9 
for C = 7 and C = 20. The experimental data represented by the symbols (taken from 
refs. 7, 10, 11, and 12) were obtained from corrugation- stiffened panels which exhibited 
flutter at Mach numbers from 1.2 to 3. The panels have a wide variety of edge-support 
conditions ranging from very weak supports, with the corrugations unsupported, to the 
attachment of the corrugations to a rigid substructure. Deflectional stiffnesses at the 
supports were not measured and were difficult to define. However, an estimate of the 
stiffness was made for the flexibility of the corrugations at the supports as well as the 
flexibility of the attachments. Details of deflectional stiffness calculations are presented 
in reference 12. 

The diamond symbols in figure 10 are for values of C between 5.5 and 7 and with 
the exception of two test points, all lie in the flutter region above the curve for C = 7. 

The square symbols represent panels with intermediate values of C between 9 and 15, 
and the circular symbols represent panels with high values of C between 19 and 24. 

The data are seen to verify the trends of the theoretical curves and illustrate the fact 
that corrugation- stiffened panels tested to date have fluttered at dynamic pressures far 
less than the value provided by the conditions for all edges simply supported. Some of 
the data appearing in figure 10, along with data on isotropic panels have been used to gen- 
erate empirical flutter envelopes in an effort to represent the lower limit of test data over 
a wide range of proportions. (See, for example, ref. 13.) No differentiation was made 
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between data for corrugation- stiffened panels (which show the large disparity from theory 
for simple support boundary conditions) and isotropic panels (which compare well with 
theory for simple support boundary conditions). It can be seen that such a flutter bound- 
ary would be overly conservative for any panel that had appreciable support stiffness. 
However, proper account of support flexibility results in reasonable agreement between 
theory and experiment. 


DESIGN RECOMMENDATIONS 

For the purpose of design of stiffened panels for flutter-free operations on the 
hypersonic vehicle, limiting flutter boundaries have been presented for both spanwise 
and streamwise orientation of maximum- stiffness corrugations. (See figs. 8 and 9.) 

Since these conservative boundaries show that support flexibility results in a rapid 
reduction in ^^r since some flexibility is inherent in the fabrication of complex 
stiffened structures, some degradation of the flutter boundary must be anticipated. 

Recently , a method of support attachment has become apparent that may provide a 
sizable increase in A cr in the presence of support deflectional flexibility. This method 
involves the addition of a local torsional stiffness along the boundaries normal to the 
maximum-stiffness corrugations. For example, a corrugated doubler strap along the 
ends of corrugations provides considerable torsional stiffness at the boundary. The 
effect of torsional stiffness has been determined by an extension of the analysis of 
con ugation-stiffened panels presented in reference 8. A sample calculation is shown 
in figure 11 where plotted as a function of Kp for Kp) = 10. For this calcu- 

lation a/b = 1, C = 7, and -A x = 100. For K-p = 0, the value of \£ r is the same as 
a point on the curve K D = 10 from figure 6. As K T increases, \£ r is seen to 
increase and the curve for Kq = 10 eventually becomes asymptotic to the curve for 
k d = 00 (that is, for all edges simply supported). Thus, the addition of torsional bound- 
ary stiffness can compensate for lack of deflectional stiffness. 

The effect of Krp on flutter has been verified to a degree by a recent experimental 
investigation conducted on corrugation- stiffened panels in the Langley 9- by 6-foot ther- 
mal structures tunnel at Mach 3. The panels were constructed with corrugated doublers 
welded to the panel corrugations at the edges. (See fig. 12(a).) These corrugated dou- 
blers provided torsional spring stiffness Kp; they were 0.60 inch (15.2 cm) wide and 
three times the thickness of the corrugated material. The panels were attached at the 
ends of each of the corrugations to carefully machined cantilevered beams (fig. 12(b)) in 
order to facilitate easy calculation of deflectional stiffness K D and rotational stiffness 
Kr. The numerical value of Kp was obtained from measured natural frequencies shown 
in figure 13. The circular symbols are the measured frequencies corresponding to one 
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half-wave in the direction of the corrugations and to the first seven half-waves in the weak 
bending direction. The dashed curve gives the theoretical frequencies for the calculated 
values, K d = 0.64 and K R = 0.15, for K T = 0. The solid curve yields theoretical fre- 
quencies for the same values of Kp and K R , but for = 1.2 and provides a realistic 
correlation of the measured frequencies, as can be seen from figure 13. 

The panel with torsionally stiffened edges was tested at Mach 3, and the flutter data 
obtained is compared in figure 14 with the theoretical boundaries from figure 10. The 
experimental data obtained from the torsionally stiffened panel is shown by the bar sym- 
bol at K d = 0.64. The panel was tested in the presence of aerodynamic heating and the 
height of the bar represents the range of flutter dynamic pressure as the result of thermal 
stress. The limited amount of data prevented definition of the experimental value of A cr 
for the unstressed panel. The theoretical prediction of A^ r (shown by the cross) is for 
zero stress and K-p = 1.2. If experimental data had been obtained at zero stress, the 
comparison would be much improved. The significant fact is, however, that in the absence 
of a doubler strap along the corrugated edges (that is, for Kt = 0), the panel would prob- 
ably flutter at a value of A£ r slightly above the theoretical curve for C = 7 (the value 
of C of the panel tested is 6.1) at K D = 0.64. The inclusion of the doubler may have 
increased A^r by as much as a factor of 30. Thus, it appears that a simple doubler 
strap at the corrugated edges or other methods of attachment which provide local tor- 
sional rigidity for the stiffening elements may result in large increases in flutter mar- 
gin and improved efficiency in panel design. 

CONCLUDING REMARKS 

Existing flutter analyses and experimental data obtained in the supersonic Mach 
number range can be applied to the design of panels for the hypersonic speed range with 
reasonable success. Examination of a hypersonic-vehicle flight trajectory indicates that 
the most critical portion from a flutter standpoint is probably in the supersonic speed 
range. 

Experimental studies in the supersonic speed range have shown that panels proposed 
for hypersonic vehicles flutter at dynamic pressures much less than anticipated based on 
classical flutter theory for panels with simply supported edges. This loss in flutter mar- 
gin is attributed to deflectional support flexibility which is inherent in the usual fabrica- 
tion of stiffened structures. Simple design curves are presented for such orthotropic 
panels and have been corroborated by the correlation of a large quantity of experimental 
data on corrugation- stiffened panels covering a wide range of test conditions and support 
conditions. Finally, it is demonstrated experimentally and analytically that even in the 
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presence of deflectional support flexibility, a large improvement in flutter margin can be 
realized by the addition of torsional restraint at the supports. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., October 9, 1969. 
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APPENDIX 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 

The International System of Units (SI) was adopted by the Eleventh General 
Conference on Weights and Measures, Paris, October 1960, in Resolution No. 12 (ref. 14). 
Conversion factors for the units used herein are given in the following table: 


Physical quantity 

U.S. Customary 
Unit 

Conversion 

factor 

(*) 

SI Unit 

Length 

in. 

0.0254 

meters (m) 

Stiffness .... 

lbf-in. 

0.113 

newton-meter (N-m) 

Pressure .... 

Ibf/ft2 

47.88 

newtons/meter2 (N/m2) 


^Multiply value given in U.S. Customary Unit by conversion factor to 
obtain equivalent value in SI Unit. 


Prefixes to indicate multiples of units are as follows: 


Prefix 

Multiple 

kilo (k) 

10 3 
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Figure lr Comparison of experiment and theory for the flutter boundary of clamped isotropic panels at various Mach numbers. 










Figure 6.- Flutter boundary for flexible supports for flow in weak direction. (Modal solution.) Numbers 


on curves designate modes. 








Figure 9.- Minimum flutter boundary for highly orthotropic panels with flow in weak direction. 








(a) Torsionaliy stiffened edge. 



(b) Beam support details. 

Figure 12.- Beam-supported corrugation-stiffened panel with torsionaliy stiffened edge. 
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Figure 13.- Comparison of experimental and theoretical natural frequencies 
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Figure 14.- Comparison of theory and experiment for flutter of torsional restrained panel. 


